Klein tunneling and Dirac potentials in trapped ions 
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We propose the quantum simulation of the Dirac equation with potentials, allowing the study of 
relativistic scattering and the Klein tunneling. This quantum relativistic effect permits a positive- 
energy Dirac particle to propagate through a repulsive potential via the population transfer to 
negative-energy components. We show how to engineer scalar, pseudoscalar, and other potentials in 
the 1 + 1 Dirac equation by manipulating two trapped ions. The Dirac spinor is represented by the 
internal states of one ion, while its position and momentum are described by those of a collective 
motional mode. The second ion is used to build the desired potentials with high spatial resolution. 



The Dirac equation describes the successful merge of 
quantum mechanics with special relativity [1] , predicting 
the electron spin and antimatter. Surprisingly, it also 
predicts some conflictive quantum relativistic effects, as 
the Zitterbewegung Q and the Klein paradox which 
have been discussed theoretically though never experi- 
mentally tested. Presently, there is growing interest in 
different aspects of quantum simulations [4, 5]. Recently, 
the quantum simulation of a free-particle Dirac equa- 
tion in trapped ions 0, Q has generated a dialogue be- 
tween relativistic quantum mechanics and quantum op- 
tics [H 0] ■ Here, we propose the quantum simulation of 
the Dirac equation with potentials in trapped ions, allow- 
ing us to explore the unintuitive physics of relativistic 
scattering, especially when compared with Schrodinger 
quantum mechanics. We also show that in 1 + 1 dimen- 
sions only the scalar and pseudoscalar potentials can be 
used to confine a Dirac particle. Moreover, we discuss 
how electromagnetic potentials give rise to Klein tunnel- 
ing (loj . in which a particle propagates throu gh a repul- 
sive potential by turning into its antiparticle 

The Dirac equation with a covariant potential [l[ can 
be written as 



—ihd — mc V r , 



rh = 0, 



(1) 



where we used Feynman's notation, JL = 7*'A M , being 
7 M the Dirac matrices and a four-vector. The exter- 
nal potential V cov can take many forms that transform 
differently under Lorentz rotations and boosts. In par- 
ticular, V cov includes scalar potentials that add to the 
relativistic mass term, as well as pseudoscalar, electric 
and magnetic potentials, among other cases. Multiply- 
ing Eq. (dJ by 7 := j3, and introducing the vector of 
matrices a = 7°7, the momentum operator p — —iKV, 
and the time t = xq /c, we obtain 



ihdtip = \ca.p + mc 2 (3 + (3V cov ] ip. 



(2) 



matrices, a and j3, determining different Dirac represen- 
tations. In the case of 1 + 1 dimensions, in particular, a 
and /3 can be chosen as any two different Pauli matrices. 
For simplicity, we will use a real representation, 



7° = = <r», Y 



(3) 



There are six ways to introduce a potential in the Dirac 
equation depending on its behavior under Lorentz 
transformations. Using the pseudoscalar operator 7 s := 
ij\ „7 M and the tensor er MI/ := 17^7", these are 



V"" + + 



(4) 



■'"*' /*„„vV'". 



The potential V transforms as a scalar and mimics an 
induced mass term. The field A transforms as a four- 
vector and corresponds to the electromagnetic poten- 
tial acting on a charge q, with the electric potential 
Aq = <fi and a three-vector component A. The remain- 
ing four potentials B, A, V and B, transform as matrices, 
pseudoscalars, pseudovectors, and pseudotensors, respec- 
tively, and behave as anomalous field moments [l|. 

In 1 + 1 dimensions the landscape simplifies consider- 
ably and we can choose a parametrization such that 



ihd tip 



~>x (p 

H(q)iP, 



<l 



A) +q(f>+ {mc 2 + V)a z - qV<7 % 



1> 
(5) 



Up to unitary transformations, there are different sets of 



where only A, <f), V, and V are nonzero. We will first 
focus on potentials that are linear in the particle position 

V = v sc x, (j) = v e ix/e, A = v mag x/e, and V = v ps x/e, 

(6) 

where e is the unit charge, and we will assume q = ±e. 

All the preceding potentials can be simulated with 
trapped ions |U . Let us consider a string of two 
trapped ions, 1 and 2, which could be 40 Ca + ions 
with long-lived internal states \Si/2,Tn = 1/2) and 
l^5/2, m = 3/2). The first ion will encode a Dirac 
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FIG. 1. Setup for simulating a Dirac particle in an external 
potential. A string of two trapped ions is manipulated with 
addressed laser beams, coupling to each ion's internal state 
and the collective motion. The Hamiltonian for a free Dirac 
particle and the potentials A, V and V can be implemented 
directly by the laser light impinging on ion 1. The electro- 
static potential <f> is implemented by red and blue sidebands 
applied to the auxiliary ion to create an interaction oc crfx, 
with the auxiliary ion prepared in an eigenstate of erf- 



spinor in those internal states, while the second ion will 
be used as an ancilla to implement potentials. Assum- 
ing the validity of the Lamb-Dicke approximation, ad- 
dressed laser beams can be used to implement interac- 
tions of the forms H c = Ml(afe 



+ a)o- e 
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cr e (carrier) 



) (red sideband) and 
,-i<t>b^ (blue sideband) on 

each ion j — 1,2. Here, fl(& ir ) and 4>(h tT ) are the Rabi 
frequency and phase of each light field, rj <C 1 is the 
Lamb-Dicke parameter and (oT) and a' (a) are rais- 
ing (lowering) operators for the ion's internal states and 
for a collective motional mode. For the one-dimensional 
case, we could use, for instance, the center-of-mass mode. 

As shown in @, 0], the Dirac Hamiltonian for a free 
particle, H { £ cc — co~ x p + mc 2 a z , can be engineered us- 
ing a simultaneous blue and red sideband interaction of 
equal strength with the appropriate relative phase, am- 



2A 



plitude and detuning Ml := mc . Here, p 
A = y^h/4rhuj is the size of the ground state wave packet, 
m is the mass of a single ion and uj is the trap frequency. 
We consider the laser coupled only to the internal state 
of ion 1 (see Fig. 1) and show how to extend this model 
to include each of the potentials V, V, A, and <j>. 

The potentials A and V can be created simultaneously 
by changing the phase and the intensity of the sidebands. 
Note that the total Hamiltonian (co~ x p — eo~ x v m agX — 
ev ps X(7y), with x = (aJ + a)A, can be decomposed in 
terms of a^a^, and aa^ . For the implementation of the 
scalar potential, we can use two additional light fields 
acting on ion 1 with Rabi frequency f2 sc , one blue de- 
tuned by lj/2 from the spinor transition, the other red 
detuned by the same amount to create a Hamiltonian 
term \hrf\} sc xa\ j uj A 13j. Alternatively, we can choose 
another representation in which the mass term acquires 



the Pauli matrix a y . In this case, the term mc 2 is created 
by a carrier excitation, whereas the scalar potential can 
be engineered from blue and red sidebands in the same 
way as A and V. Finally, the electrostatic field can be 
implemented by driving a blue and red sideband simul- 
taneously on the auxiliary ion. With appropriate phases 
the resulting Hamiltonian term becomes cx erf x. Prepar- 
ing the auxiliary ion in an eigenstate of erf, this operator 
can be replaced by its eigenvalue and the interaction re- 
duces to the desired form. A general Hamiltonian for this 
system that can be compared to Eq. ([5]), in the rotating 
frame and after rotating-wave approximation, reads 

H = hr) (n b e 1<t ' b a t a+ + n r e^ r aa^ + H.c. 

+h V n 2 a^x/A + h(n + 4r]n 2 sc x/ujA)al (7) 

Here, the first two terms with their Hermitian conjugates 
describe blue and red motional sidebands, the term in- 
volving the second ion describe a conditional displace- 
ment, and the last term involves constant and position- 
dependent Stark shifts. The relations between the Dirac 
model in Eq. ([5]) and the ion system in Eq. ([7]) are 



m = 



hr]tt 2 /A = v e i, Ah-qVt 2 sc /uA 



(8) 



while the relative weights of c, v mag , and v ps , can be set 
by changing phases 4>b,r and amplitudes fl r ,b- 

The creation of an electrostatic potential via a detuned 
laser acting on ion 2 yields the interesting case of inter- 
action terms hij^xa^ / A + h^a^- For a large detuning 
0,2 r]Cl2 this interaction becomes effectively cx i 2 crf , 
allowing the simulation of quadratic potentials [IH • 

Electric potential and Klein tunneling.- Not all po- 
tentials can actually confine a Dirac particle. To prove 
this let us introduce an antiunitary operation known as 
charge conjugation ip c = Kip = Cip* , combining com- 
plex conjugation and a unitary matrix C7 M *C _1 = — 7 M . 
The charge conjugate spinor satisfies a Dirac equation 
with opposite charge ihdtijjc = H(—q)tp c , in which cer- 
tain terms of the covariant potential which are propor- 
tional to the charge "q" have changed sign. We can also 
show that the negative energy branch of a Dirac equation 
can be antiunitarily related to the positive energy branch 
with opposite charge KH(q)K^ 1 = —H(—q). In other 
words, the negative energy branch of the Dirac spectrum 
are indistinguishable from particles with opposite charge. 
Moreover, these antiparticles will see certain components 
of the covariant potential, A, V, A, with opposite sign as 
for their positive energy counterparts. 

The fact that positive and negative energy states see 
different effective potentials, and that both components 
are coupled, allows the apparition of the so-called Klein 
paradox The electric potential in Eq. ([6]) is repul- 
sive for positive-energy charged particles in a large re- 
gion, (q/e)v e ix > 0, but it will allow the antiparticle or 
negative energy states tunnel into a region that would be 
forbidden in the (Schrodinger) nonrelativistic regime (see 
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FIG. 2. Klein tunneling in a repulsive potential, <f> e i(x) = 
v e ix. (a) In position space, the particle may bounce back (I) 
or enter the forbidden region by reducing its kinetic energy 
(II). (b) In momentum space, (I) and (II) correspond to the 
particle having positive energy or turning into an antiparticle. 

Fig. [2]). The situation is even more interesting, for if wc 
send a massive charged positive-energy particle against 
an electrostatic potential barrier and the particle has 
enough energy, it will split into a positive and negative 
wavepackets, the former bouncing back and the latter 
penetrating into the energy barrier (see Fig. [3]). 

We explain now how Klein tunneling can be interpreted 
in terms of a quantum optical concept: the Landau-Zener 
tunneling [la ]. Working in momentum space, where x = 
(+ihd p ), and with charge q — e, the electric potential 
v e ix is equivalent to a deceleration or a decrease in the 
particle momentum, which can be compensated using the 
change of variables ip(p, t) = £(p + v e it), 

ihd t £, = [ca x (p - v e it) + mc 2 a z ]^. (9) 

This equation corresponds exactly to a Landau-Zcncr 
process in which an effective magnetic field along the x 
direction is increased linearly in time. 

The dynamics is summarized in Fig. [2j If the particle 
initially moves against the slope, that is p > 0, the poten- 
tial will decelerate the particle which may (I) bounce back 
or (II) tunel inside the potential barrier. The first case 
corresponds to adiabatic transition in the Landau-Zcncr 
picture, it happens when the change of the momentum 
is small, fic|t> e ;| <C m 2 c , and is the typical behavior of 
nonrclativistic particles. However, for large enough slope 
or in the relativistic limit, it becomes possible for the 
particle to switch branch, acquiring negative kinetic en- 
ergy and entering the originally forbidden region. More 
precisely, the probability of transition can be computed 
with the Landau-Zener formula as follows 

P " = CXP H^J- (10) 

Note that once the particle has switched to negative en- 
ergy states, it has an opposite charge and the potential 
4> e i accelerates the particle indefinitely in the opposite 



direction. This makes such a potential effectively a filter 
that separates particles from antiparticles. 

Using Eq. ||!JJ|, it is easy to perform numerically accu- 
rate simulations of the Klein scattering for arbitrary ini- 
tial conditions and different ratios of the potential slope, 
m 2 c 4 /hcvei. In Fig. [31 we see that for very small masses 
the particle is either not confined, or it splits into a nega- 
tive charge and a positive charge components [Fig.[3jD] . In 
the nonrelativistic limit, m 2 c A ^> hcv e i, though, the par- 
ticle mostly bounces back from the barrier. The first plot 
corresponds only to situation II, the second plot combines 
I and II, and finally the third situation is just I. 

Scalar potential- As we discussed above, the Klein 
paradox only applies to potentials which are not invari- 
ant under charge conjugation. A scalar potential fisj ] 
such as V in Eq. (Q}, acting similarly on both particle 
and antiparticle states, can still confine and reproduce 
the physics we are used to in the nonrelativistic limit. 
In particular for our choice of linearly growing potential 
V = v sc x, with v e i,v magi v vs = 0, it is possible to show 
the existence of bound orbits. This is best analyzed by 
squaring the effective Hamiltonian, 

H 2 = cp 2 + (v sc x + mc 2 ) 2 + chv sc o-y. (11) 

The first implication is that the scalar potential eigenen- 
ergies correspond to those of a harmonic oscillator 

E 2 = 2hcv sc + i ± l^j . (12) 

The second implication is that the conserved quantity H 2 
is defining elliptical orbits in the phase space of position 
and momentum once more as a harmonic oscilla- 

tor but these orbits are now centered around the point 
x = mc 2 /v sc . These closed and bound orbits exist both 
for the Dirac particles and antiparticles, implying that 
both are confined around the same trajectories. 

Pseudoscalar potential and Dirac oscillator.- The last 
analyzed case is the pseudo-scalar term v ps xa y with 
v ps = mujc. There is a similarity between 

H = c(a x p + (T y mujx) + mc 2 a z (13) 

and the Hamiltonian of a 1 + 1 Dirac oscillator |8|, 

ihdt^P = \ca{p + imuifix) + mc 2 0\ (14) 

but this coincidence only happens because of our re- 
stricted dimensionality, 7 5 = a(3. With the oscillator 
length scale, a osc = yjh/mui, and the Rabi frequency 
Ml = V Humid 2 , the pseudoscalar potential can be refor- 
mulated as a detuned Jaynes-Cummings model 

ihdt^p — [hfl(ia — ia + a) + mc 2 o~ z ~\ ip, (15) 

which should be easy to simulate. A similar calculation 
as before produces again closed phase space orbits 

p 2 + x 2 mu) + ha z = E 2 , (16) 

mui 
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which are a consequence of the JC discrete spectrum 



E n = ±mc 2 \ n 7t + 1 ~ mc 2 + —ntoh. (17) 

V mcr 2 

Note also that a replacement of ui — > —u is still a Dirac 
oscillator, but then it turns into an anti-JC Hamiltonian. 

Conclusions and outlook.- We have proposed a quan- 
tum simulation of the 1 + 1 Dirac equation with poten- 
tials in trapped ions, bringing together the physics of rel- 
ativistic quantum mechanics to a controllable tabletop 
experiment in quantum optics. Our simulation protocols 
can be combined with tools for monitoring and preparing 
the quantum state of the ions. In particular, newly de- 
veloped techniques 0, E3] make it possible to obtain the 
position, momentum and probability distribution |^(x)| 2 
of the particle in an efficient way, allowing the "frame- 
by-frame" reconstruction of the scattering event. Besides 
this, the initial state of the Dirac particle can be accu- 
rately engineered in position, momentum and even en- 
ergy branch This is important for the electrostatic 
potential because, by preparing a particle with positive 
energy, we could see a full reflection for shallow slopes 
and full transmission for steeper slopes, making the Klein 
paradox visible. As an example, similar to Ref. [HI, one 
could prepare a positive energy wave packet with average 
momentum p — 4ft/ A and choose experimentally accessi- 
ble parameters £1/, = fi r = 2ttx 20 kHz, Qi = 2ttx 1 kHz, 
i] = 0.05 and 0,2 = 2-7TX 50 kHz. This permits to obtain 
a tunneling probability of 0.5 which could be observed in 
the lab. Numerical simulations show that the whole tun- 
neling dynamics with these parameters takes place within 
1 ms, which is well within the motional and internal state 
coherence time. To study the Klein tunneling, see the 
case of electric potential and Fig. [3l the slope strength 
of the mass of the simulated particle can be varied by 
changing laser intensities and frequencies. 
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